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We explore theoretically impurity states in the antiferromagnetic spin-density wave state of the 
iron arsenide. Two types of impurity models are employed: one has only the intraband scattering 
while the other has both the intraband and interband scattering with the equal strength. Inter- 
estingly, the impurity bound state is revealed around the impurity site in the energy gap for both 
models. However, the impurity state is doubly degenerate with respect to spin for the first case; 
while the single impurity state is observed in either the spin-up or spin-down channel for the second 
one. The impurity-induced variations of the local density of states are also examined. 

PACS numbers: 71.55.-i,75.30.Fv,75.10.Lp 



The recent discovery of iron-based superconductors [l| 
has triggered intensive efforts to unveil the nature of and 
interplay between magnetism and superconductivity in 
this family of materials. Series of iron arsenide have 
been synthesized, vifhich possess many similar features 
of the normal and superconducting states. Experimental 
measurements have reported that the undoped ReFeAsO 
(where Re= rare-earth metals) and AFe2As2 (where 
A=divalent metals such as Ba, Ca, Sr) compounds ex- 
hibit a long-range antiferromagnetic spin-density-wave 
(SDW) order 0, i, i, H, H, 0- Upon electron/hole dop- 
ing the SDW phase is suppressed and superconductivity 
emerges with up to above 50 K [1, i, [H El, [13 ■ 

At present, there is likely certain controversy on 
the understanding of the SDW state of the undoped 
FeAs-based parent compounds. Two kinds of theories 
have been put forward: 1) the itinerant antiferromag- 
netism, which takes advantage of proper Fermi surface 
(FS) nesting (or strong scattering) between different FS 
sheets IJ, [iJ, |l5|, ll6|; and 2) the frustrated Heisen- 
berg exchange model of coupled magnetic moments of 
the localized d-orbital electrons around the Fe atoms 
17 , 18 . [Til [ioj . As for the itinerant electronic behav- 
ior, first principle band structure calculations [2l| based 
on the density functional theory (DFT) indicate up to 
five small Fermi pockets with three hole-like pockets cen- 
tered around the F point and two electron-like ones cen- 
tered around the M point of the folded Brillouin zone 
of the FeAs layers, which have partially supported by 
the angle-resolved photoemission spectroscopy (ARPES) 
from different groups [IS, S [13, S, IH- Motivated by 
the DFT calculation and experimental measurements, in 
Refs. jl5l. the excitonic mechanism 27 1 of itinerant 
carriers are employed taking account of the FS nesting 
between electron and hole pockets and the SDW phase 
are associated with triplet excitonic state, which can be 
understood as condensate of triplet electron-hole pairs 



27| 



In this paper, we explore theoretically the effect of a 
single impurity on the local electronic structure of an Fe- 
based antiferromagnet in the triplet excitonic phase. It 
is shown that impurity bound states are formed inside 
the SDW gap, which may be observed experimentally by 
local probes. Before introducing the impurity, we first 
propose an effective model Hamiltonian to address the 
triplet excitonic state, 
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where T — (0,0), X — (7r,0). We use the index i to 
label different valence bands around F point. Around 
X and Y points, there are two conduction bands, diko- 
and Ck+Xo- are the annihilation operators of electrons in 
the Ti and X bands. Theoretically X and Y are two 
equivalent nesting directions. Note that, the structural 
phase transition occurred just above/on the SDW transi- 
tion breaks this equivalency. Without loss of generality, 
it is assumed that only conduction band around the X 
point couples with the valence bands around the F point, 
which is characterized by the mean-field order parame- 
ters AiCTo-'- For the triplet excitonic phase (SDW), we 
have real order parameters satisfying = ^^iii a-i^d 
A,ti = A^T = [13. 

eri(k) and ex(k) are used to denote the band disper- 
sions of the nonmagnetic normal state. For k in the vicin- 
ity of the F point (therefore, k -I- X in the vicinity of the 
X point), the normal-state energy dispersions have ap- 
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FIG. 1: Schematic plot of (a) the Fermi surfaces; and (b) 
the band dispersions of the valence (hole) band and conduc- 
tance (electron) bands in the unfolded Brillouin Zone for the 
undoped parent compound. See text for detail. 



For the reason that the order parameter Ai is set to 
zero, there is no couphng between the Fl band and the 
X band. The Hamiltonian of Eq. ([1]) is reduced to a 
model of two bands with one valence band (F2 band) 
and one conduction band (X band). Introducing the two- 
component Nambu operator, ■0^^.^ = (rf2kcr' "^Ic-i-Xct)' 
model Hamiltonian can be simplified as 
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where an impurity term has been added with the form, 

iiimv = X! V^L^k,k'V'k'<T, (8) 



k,k' 



proximately the 2D parabolic forms 



er.(k) : 
£x(k + X) = 



2mri 



^0 ' 
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as schematically shown in Fig. [TJ Here mn and mx 
are the corresponding effective masses. In describing the 
X band, the elliptic FS is approximated by the circu- 
lar one for simplicity, eg' (ej^) denotes the top (bot- 
tom) of the hole (electron) bands. According to the 
ARPES measurement 122] , two hole-like Fermi pockets 
are revealed around the F point for undoped BaFe2As2. 
The band parameters extracted from the experimental 
data are as follows, mri ~ 2.8me, mr2 ~ lAnie, 
and mx ~ 6.5me, where rUe is the mass of bare elec- 



tron. 



.n 



4 meV, eP w 16 meV, and £{f 



24 



meV. These parameters indicate that the nesting be- 
tween the F2 band and X band is much better than 
that of the Fl band. Therefore it is natural to assume a 
larger order parameter A2 and a vanishingly small Ai. 



Let Eg 
E, 



(eo + Co )/2 and mo = (' 



T2 




e^)/2. Here 



G = —2Eg denotes the indirect gap between the top of 
the F band and the bottom of the X bands. Therefore, 
Eg > describes a semimetal and Eg < a semiconduc- 
tor. With the help of Eg and ^0 and a further assumption 
of mr2 = "mx = m « 7me. we can re-express the energy 
dispersions as 



£r2(k) = -e(k) - ^0, 
£x(k + X) = £(k)-^o, 

£(k) = |ik2-i?^ 

2m 
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(4) 
(5) 

(6) 



Note that for ^0 = 0, the hole and electron bands are per- 
fectly nested since £r2(k) — — ex(k -I- X) and the system 
is unstable with respect to infinitesimal Coulomb inter- 
action while for nonzero ^0 finite strength of Coulomb 
repulsion is needed. 



where ?7k,k' represents a 2 x 2 matrix of the scattering 
potential associated with non-magnetic impurities. Here 
we use Ao- to denote A„„ for short. The Green's function 
method is applied to study the single impurity effect. The 
matrix Greens functions are defined as 



C(k,T;k',r') = 

G'crcr(k,k',iLJ„) = 



(T.[^k.(r)4'.M)), (9) 
dTG,,(k,r;k',0)e"^"^ (10) 



uj + iO+). (11) 



From the Hamiltonian defined in Eq. ^ we can derive 
the bare Green's function 



w-er(k) -Ao- 
-A^ tj-ex(k) 

Qtq + A0.f1 - £(k)-r3 



(12) 
(13) 



where lu = uj + fiQ. tq is the 2x2 unit matrix, and fi^a are 
the pauli matrices. The T-matrix approximation is em- 
ployed to compute the Green's function in the presence 
of impurities. For a single impurity, the T-matrix exactly 
accounts for the multiple scattering off the impurity. The 
single-particle Green's function G can be obtained from 
the following Dyson's equation, 

G„,(k,k',c.) = G^,(k,c.)Vk'+G°,(k,c.) 

f,,(k,k',c.)G°,(k',c.), (14) 

where the T matrix is given by 

f,,{k,-k',uj) = t/k,k'+5I'^k,k"GL(k",^)ro.(k",k',L^). 

k" 



(15) 

For a point-like scattering potential interacting with itin- 
erant carriers just on the impurity site, the scattering 
matrix is isotropic, C/k,k" = U. The above equation is 
greatly simplified 



(16) 
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where G^^{uj) 
we obtain 



where 



J2k ^aa 0^^ ^) ■ After some derivation 



= {ujro + A^fi) + 7(u))t3 



(17) 



with f' = To + Ti. The energy of the impurity bound 
state is again determined by the pole of Ta-a-{uj). From 



det[2V~lj, - T'G°„{n)] = we have the equation for fi, 



n + A. 



(24) 
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with Ec denoting the high-energy cutoff and A^o — 
ma? / the density of states per band per spin. Note 
that a{Cj) and 7(0)) are independent of the spin index a. 

The first impurity model we study is the scattering- 
potential matrix with only intraband scattering terms, 
i.e. ij — VimpTo, which was adopted in Ref. [28[ to study 
effect of many impurities. From Eq. (|16p . we obtain 

TM - [ro-UG'M]-'U = [VrJ^-GU^)]-\ (20) 

The energy of the impurity bound state is determined by 
the pole of f^^iuj), determined by det[V-^l - G°<,(0)] = 
0. Setting c = (TrA^o^mp)"^? we have the equation for 
the energy of impurity bound state 

" 2c-i£l£= - a(f^)2 - 7(f])2 = 0. (21) 

Ai-n-^ 



For the spin triplet excitonic phase A| = — A|, the 
above equation gives rise to impurity states with the same 
bound energy, i.e. the impurity states are doubly degen- 
erate. Generally, the above equation has to be solved nu- 
merically to obtain the bound energy Q. However, we can 
get some analytic results under certain approximations. 
Under the wide-band approximation Ec,Eg ^ |Act|, 
q;(0) « 1 and "f{^) » 7o = t^^^ ^''^{Ec/Eg), so we have 



ia: 



sgn(,cj 



2 
75 



V(l - C2 + 7^)2 + 4c2 ' 



(22) 



and furthermore if the system has approximately the 



particle-hole symmetry Ec 



Eg, then 70 w and 



n/\Ac\ = sgn(c)(l — c2)/(l-|-c2) from the above equation. 

For the second impurity model, the four matrix ele- 
ments of U is assumed to be the same, i.e. the intra- 
and inter-band scattering terms are the same [2^ with 
U = VimpiTQ +fi)/2. Then the T-matrix according to 
Eq. (dll) is 



r..(^) = [2v; 



1 

imp 



(23) 



From the above equation we find that J7 is independent 
of the function of 7(0)) for this case, which reflects the 
particle-hole asymmetry. Before solving the above equa- 
tion for the bound energy, we study the existence of the 
impurity state. Because < A2 , Cl + has the same 
sign as that of A^.. Therefore, the solution of Eq. 
exists only if the sign of c is opposite to that of A^r- For 
the SDW state, i.e. the triplet excitonic phase, we have 
A| = — A| and so there is exactly one impurity bound 
state in either the spin-up or spin-down channel. We 
may assume A^ = — A| = A > as well, then for attrac- 
tive scattering Vimp < 0, the impurity bound state only 
exists in the spin-up channel and its energy is given by 
Q/A= — (1 — c2)/(l-fc2) under the wide-band approxi- 
mation. If Vimp > 0, however, the impurity state will be 
in the spin-down channel, and Sl/A = (1 — c2)/(l -|- c^). 
In general, the impurity bound-state energy is given by 
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a: 



sgn(c) 
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(25) 



in the valid regime of the wide-band approximation. 

To apply the theoretical results to the iron arsenide, 
we try to pin down the parameters of our model by ex- 
tracting them from the available experimental data for 
BaFeaAsa ;225. 



that En 



20 meV. 

-1 



« 16 meV and €q 
mr ~ mx ~ 7.0 me 

At = -A4 



24 meV so 
and there- 
fore No w 1.2 eY-\ A| = -Aj, = A « 20 meV. The 
high-energy cutoff is set as Ec = 500 meV, which is of the 
same order of magnitude as the band width. Note that 
Eg extracted from experimental data is very small, which 
is in the same order of magnitude of the order parameter 
A. Therefore, neither the wide-band approximation nor 
the particle-hole symmetry can be applied to the present 
case. Eqs. (l^lT) and have to be numerically solved. 

Now we examine the local characteristics induced by 
the impurity by looking into the variation of the local 
density of states (LDOS), which can be probed by the 
scanning tunneling microscopy (STM). The LDOS is de- 
fined as 



= -- V/77i{rr[G,,(r,r',c.)]}, 



(26) 



where Go-CT(r, r', t^) the Green's function in real space. 
Applying the T-matrix approximation we have. 



G,,(r,r',c.) - G^,(r,r',c.) + GO,(r,0,L.) 
f(c.)GO,(0,r',c.), 



(27) 
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Substituting Eq. (|27|) into Eq. ([26)) we may single out the 
variation of LDOS due to the presence of the impurity 
potential, 



GL(0,r,c.)]}. (28) 

For the second impurity model,30], Fig. [2^ a) shows the 
LDOS as a function of energy Co on the impurity site, 
namely iV(0,a)), while Fig. E^b) the impurity-induced 
LDOS at the bound energy as a function of radial dis- 
tance r off the impurity site, i.e. Ni„ip{r,n). Vimp has 
been set as -0.36, -0.6, and -4.0 eV, giving rise to the 
impurity bound states seen as the sharp peaks located 
respectively at the energies f2/A = 0, —0.5, and —0.99 in 
Fig.HKa). The probability densities of these bound states 
exhibit a kind of exponential decay with the Friedel os- 
cillation, as seen in Fig. [Ifb). Introducing two length 
scales, ^1 and ^2 to characterize the oscillation and de- 
cay, we obtain the asymptotic behavior of Nimp{r, ^) for 
large r, 



Nimp{r,^l) ocr ^ cos^(7rr/^i) exp(-r/^2)- 



(29) 




FIG. 2: LDOS on the impurity site as a function of energy uj, 
i.e. N{0,il}) (a), and impurity-induced LDOS at the bound 
energy f2 as a function of radial distance r off the impurity 
site, i.e. 7V(r, SI) (b). r is in unit of a with a the lattice 
constant of Fe-Fe plane. Black, red, and blue lines correspond 
to Vimp = —0.36, —0.6, —4 eV, respectively. 
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